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Abstract 

Phenotypic plasticity and genetic adaptation are predicted to mitigate some of 
the negative biotic consequences of climate change. Here, we evaluate evidence 
for plastic and evolutionary responses to climate variation in amphibians and 
reptiles via a literature review and meta-analysis. We included studies that either 
document phenotypic changes through time or space. Plasticity had a clear and 
ubiquitous role in promoting phenotypic changes in response to climate varia- 
tion. For adaptive evolution, we found no direct evidence for evolution of 
amphibians or reptiles in response to climate change over time. However, we 
found many studies that documented adaptive responses to climate along spatial 
gradients. Plasticity provided a mixture of adaptive and maladaptive responses to 
climate change, highlighting that plasticity frequently, but not always, could ame- 
liorate climate change. Based on our review, we advocate for more experiments 
that survey genetic changes through time in response to climate change. Overall, 
plastic and genetic variation in amphibians and reptiles could buffer some of the 
formidable threats from climate change, but large uncertainties remain owing to 
limited data. 



Introduction 

Climate change is forcing species to pursue their climatic 
niche as it shifts across landscapes (Chen et al. 2011). Some 
species will not track their optimal climates because they 
disperse poorly or they cannot disperse across natural or 
human barriers (Loarie et al. 2009; Schloss et al. 2012; 
Urban et al. 2013). These species face the threat of high 
extinction rates (Thomas et al. 2004; Malcolm et al. 2006; 
Schloss et al. 2012). However, phenotypic plasticity and 
genetic adaptation could allow these species to persist in 
their current locations despite climate change. Here, we 
review studies that examine evidence for plastic or genetic 
changes in amphibians and reptiles in response to climate 
variation. 

Amphibians and reptiles might be particularly sensitive 
to climate change. Both taxa have already experienced 
extensive declines and extinctions worldwide (Gibbons 
et al. 2000; Houlahan et al. 2000; Blaustein and Kiesecker 
2002; Sinervo et al. 2010). Both groups have been called 
'canaries in the coal mine' based on their supposed sensi- 
tivity to environmental change (Kerby et al. 2010; Mitchell 



and Janzen 2010). A global assessment further suggests that 
amphibians are declining faster than birds or mammals 
(Stuart et al. 2004). Concurrently, Huey et al. (2010) call 
lizards 'the new amphibians' because of their elevated 
extinction risk (Gibbons et al. 2000; Sinervo et al. 2010). 
Climate change has contributed to these threats (Houlahan 
et al. 2000; Carey and Alexander 2003; Stuart et al. 2004), 
and correlative climate envelope models predict that cli- 
mate change will cause extinction in 12-47% of endemic 
frogs and 11-49% of endemic reptiles (Thomas et al. 
2004). Sinervo et al. (2010) predicted that 20% of lizard 
species will become extinct by 2080, assuming a limited 
potential for evolutionary responses. Overall, both amphib- 
ians and reptiles are declining, and more tenuous analyses 
link their declines to climate change. Climate change will 
likely decrease fitness substantially in some populations 
and thus could generate strong selection on climate-related 
traits. However, considerations of these threats rarely 
account for the mitigating effects of phenotypic plasticity 
and evolutionary adaptations. 

Adaptive evolution and plasticity offer two ways species 
can sustain fitness despite climate change (Holt 1990; 
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Urban et al. 2012). Biologists continue to debate the 
importance of these mechanisms in mediating climate 
change impacts (Gienapp et al. 2008; Hoffmann and Sgro 
2011). Merila and Hendry (2014) argue in this issue that 
neither adaptive evolution nor plasticity should be consid- 
ered the null hypothesis for climate-related phenotypic 
change. We agree and further suggest that we should condi- 
tion our predictions about their relative importance based 
on available empirical knowledge. Toward this end, we 
evaluate the degree to which amphibians and reptiles 
respond to climate variation through evolved genetic differ- 
ences and environment-induced plasticity. We review stud- 
ies that evaluate evolutionary or plastic response to spatial 
or temporal variation in climate and perform a meta-analy- 
sis on the proportion of traits showing climate-related dif- 
ferences. We also assess if those changes are adaptive and 
how accurately phenotypic differences map to climate 
change versus other correlated environmental variables. 

Adaptive evolution and plasticity in amphibians and 
reptiles 

We now briefly outline the current understanding of evolu- 
tion and plasticity with respect to climate responses in rep- 
tiles and amphibians before conducting our literature 
review and analysis. Plasticity occurs often in amphibians 
and reptiles - to such a degree that we expect plasticity in 
almost every trait we measure, ranging from behavior to 
morphology to life-history traits (Relyea 2001; Booth 2006; 
Urban 2010). Phenological shifts are especially common in 
amphibians (Beebee 1995; Gibbs and Briesch 2001; Todd 
et al. 2011), where changes in breeding phenology are two 
to four times stronger than responses in other taxonomic 
groups (Parmesan 2007). For instance, the dwarf salaman- 
der Eurycea quadridigitata migrates to breeding ponds 
76 days later than it did 25 years ago (Todd et al. 2011). 
These phenological changes likely involve plasticity. How- 
ever, the evolutionary contribution has never been directly 
measured, making it difficult for one to conclude about its 
importance. Lizards are expected to respond to climate 
change by altering sex ratios, habitat choice, and hatchling 
traits (Booth 2006). Scientists have argued both for (Gvoz- 
dik 2012) and against (Huey et al. 2009) the importance of 
plastic acclimation to warmer temperatures in buffering cli- 
mate change effects in lizards. Huey et al. (2009) suggested 
that few tropical forest ectotherms demonstrate acclima- 
tion. In contrast, Gvozdik (2012) argued that plastic accli- 
mation to higher temperatures might be common, but 
experimental approaches used to estimate plasticity often 
are conducted over extremely limited time spans or ignore 
additional sources of plasticity. 

Can relatively long-lived amphibians and reptiles geneti- 
cally adapt to climate change? One study parameterized a 



model of evolutionary responses to climate change based 
on estimates of phenotypic variation in critical thermal 
maxima and average evolutionary rates in amphibians and 
reptiles (Skelly et al. 2007). The authors predicted that crit- 
ical thermal maxima could evolve 3.2°C in 50 years - fast 
enough to track climate changes. However, this estimate 
ignores potential limits to additive genetic variation at 
extreme temperatures or selection that is so strong that 
populations become extirpated before adaptation. 

In contrast to amphibians, scientists have argued that 
reptiles are less likely to evolve temperature-related traits. 
This more pessimistic view stems from reptiles' longer gen- 
eration times relative to rapid climate change, low herit- 
abilities, and genetic trade-offs among related traits (Janzen 
1994; Schwanz and Janzen 2008; Sinervo et al. 2010). 
Sinervo et al. (2010) estimated that recently extirpated liz- 
ard populations had faced a standardized selection intensity 
of 0.34 on average. This strong selection intensity might be 
sustained under controlled conditions in a laboratory, but 
might often lead to extirpation in wild populations. How- 
ever, empirical evidence remains limited in scope. For 
example, the evidence for low heritability in traits related 
to lizard persistence comes from a single population evalu- 
ated in the laboratory (Sinervo et al. 2010), and yet we 
know that heritabilities vary across species, populations, 
and environments. Also, genetic variation must be high 
enough to allow adaptation to rapid climate change. 
Genetic variation exists within and among populations for 
temperature-dependent sex determination in reptiles, but 
modeling studies suggest that it might not be sufficient to 
counter predicted levels of climate change (Morjan 2003; 
Mitchell and Janzen 2010). Hence, previous work suggests 
divergent views on the ability of amphibians and reptiles to 
adapt in response to climate change, setting the stage for a 
more synthetic analysis of the literature. 

Literature review and meta-analysis 

To fill the gap in our understanding of evolutionary and 
plastic responses to climate variation, we searched Web of 
Science for articles using the following search terms: 
'amphibia*/reptile*' AND ('climate' OR 'global warming') 
AND ('plastic*' OR 'evol*' OR 'adapt*' OR 'selection' OR 
'reaction norm' OR 'genotype by environment' OR 'GxE' 
OR 'phenotyp*'). We further searched literature citations 
of these articles for additional references. We read each 
article and only considered studies that presented empirical 
evidence for phenotypic or genotypic differences associated 
with climate change (e.g., temperature, precipitation). For 
each study, we recorded the climate change factor, pheno- 
typic response (e.g., laying date, size, survival), and evalu- 
ated evidence for evolved or plastic differences for each 
species or population. We found 72 estimates recorded in 
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30 studies that estimated climate-related trait variation in 
50 amphibian species and 29 estimates recorded in 24 stud- 
ies conducted on 18 reptile species (Table 1). 

We considered studies that evaluated phenotypic or 
genotypic differences in response to temporal (allochronic 
studies) and spatial variation (synchronic studies) in cli- 
mate. Allochronic studies evaluate phenotypic responses to 
climate change over time and provide the strongest infer- 
ences (Hendry and Kinnison 1999; Merila and Hendry 
2014). The alternative to the allochronic study is the syn- 
chronic one, which documents phenotypic variation in 
space and requires that space can be substituted for time. 
The space-for-time substitution offers weaker inferences 
especially for genetic change. For instance, a population 
that adapted to climate variation in the past might not 
adapt to future climate change if genetic variation is lack- 
ing. Substituting space for time likely provides more robust 
predictions for phenotypic plasticity; plasticity that is 
induced by spatial climate variation is also likely to be 
induced by climate variation through time. 

We only uncovered one study of allochronic genetic 
change in response to climate change in amphibians or 
reptiles and this one did so indirectly (Phillimore et al. 
2010). Another study evaluated genetic studies through 
time, but not necessarily on the same populations (Gibbs 
and Karraker 2006). Two other studies evaluated spatial 
differences among populations that had experienced diver- 
gence in temperature over time and thus represents a 
more robust case for applying a space-for-time substitu- 
tion (Skelly and Freidenburg 2000; Freidenburg and Skelly 
2004). Sixty- three percent of amphibian estimates and 
38% of reptile estimates were based on variation across 
space rather than time. To understand if these disparate 
designs affected overall outcomes, we tested for an effect 
of study design. 

We performed a meta-analysis on the proportion of 
traits that contributed significantly to genetic adaptation, 
plasticity, and genotype-by-environment interactions. We 
logit-transformed proportions and weighted each study by 
the inverse expected variance or np/(l-p) where n is the 
number of traits analyzed in each study and p is the pro- 
portion of traits showing a significant difference (Lipsey 
and Wilson 2000). We added 0.5 to studies with zero pro- 
portions in order to include them (Veichtbauer 2010). We 
counted the responses of different species and different 
populations reported in the same study as separate esti- 
mates. However, we included study, species, and popula- 
tions as random effects to control for their 
nonindependence in analyses. 

We conducted analyses in R (v. 3.0.0) and used the lmer 
function to conduct mixed-effects models. The significance 
of random effects was assessed using a likelihood ratio test 
(Pinheiro and Bates 2000). We also tested for effects of 



study design (allochronic versus synchronic) and used the 
pvals function in R to calculate significance values. 

Traits 

We used a word cloud algorithm (www.wordle.net) to 
visualize the frequency with which researchers measured 
different traits. This algorithm sets word frequencies pro- 
portional to font sizes. We arranged traits alphabetically 
from left to right. We included all traits analyzed across all 
studies while also consolidating traits into more general 
categories (e.g., 'breeding date' represents peak calling date, 
first date of arrival, and observation of egg deposition). 

Most studies on amphibians focused on phenological 
changes in breeding (Table 1; Fig. 1A; 31% of traits ana- 
lyzed). Research also focused on amphibian life-history 
traits related to larval growth (14%), survival (9%), and 
embryonic development (8%). We found fewer measures 
of direct physiological effects of temperature (Skelly and 
Freidenburg 2000), resistance to climate-related diseases 
(Chatfield and Richards-Zawacki 2011; Forrest and Sch- 
laepfer 201 1), or adult traits other than breeding phenology 
(Tryjanowski et al. 2006; Reading 2007; Alho et al. 2010), 
despite the potential importance of these traits in mediat- 
ing climate threats (Carey and Alexander 2003; Corn 2005; 
Pounds et al. 2006). Reptile studies focused on a wider 
range of traits including adult and offspring size (19%), sex 
determination threshold and sex ratios (15%), nesting site 
choice and depth (9%), incubation period (8%), and sur- 
vival (8%; Fig. IB). 

Evidence for plastic responses to climate variation 

The inverse-variance weighted proportion of plastic 
responses to climate variation in amphibians was 0.71 
(95% confidence interval: 0.60, 0.80; Fig. 2A) based on 69 
estimates. The random effects of species, study, and popu- 
lation did not significantly explain additional variation. 
The unweighted model indicated a similar proportion 
(= 0.70). Study design (allochronic versus synchronic) did 
not significantly alter the estimated proportion of plastic 
responses to climate variation (P = 0.220). 

The inverse-variance weighted proportion of plastic 
responses to climate variation in reptiles was 0.69 (0.53, 
0.81; Fig. 2A), based on 26 estimates. We detected a signifi- 
cant random effect of species (Zdf=i = 6.2, P = 0.012) but 
not of study or population. The unweighted model indi- 
cated a similar proportion (= 0.70). Whether the study 
design was allochronic or synchronic did not significantly 
alter the estimated proportion of plastic responses to 
climate variation (P = 0.722). 

Many of the collected studies evaluated how amphibians 
and reptiles adjusted their breeding phenology in response 
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Table 1. Summary of 26 studies on reptiles and 32 on amphibians designed to examine plastic and genetic responses of traits driven by climate 
variation. 

Common Name Species Trait Type Factor Genetic Plastic GxE Adapt Cause Time Reference 
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Common Name 


Species 


Trait Type 


Factor 


Genetic 


Plastic GxE 


Adapt 


Cause 


Time Reference 
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Common Name 


Species 


Trait Type 


Factor 


Genetic 


Plastic 


GxE 


Adapt 


Cause 


Time 


Reference 
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Critical thermal 




Y(2,5) 






Y(5) 


Y(1,2) 




(2000) 






maximum 


















Western toad 


Bufo boreas 


Breeding date 


T 




N(4) 








FD 


Blaustein et al. 


Western toad 


Bufo boreas 


Breeding date 






Y(4) 




Y(5) 


Y(1,2) 


FD 


(2001) 


Western toad 


Bufo boreas 


Breeding date 






Y(4) 




Y(5) 


Y(1,2) 


FD 




Cascades frog 


Rana cascadae 


Breeding date 






N(4) 








FD 




Cascades frog 


Rana cascadae 


Breeding date 






Y(4) 




Y(5) 


Y(1,2) 


FD 




Fowler's Toad 


Bufo fowleri 


Breeding date 






N(4) 








FD 




Spring peeper 


Pseudacris 


Breeding date 






Y(4) 




Y(5) 


Y(1,2) 


FD 






crucifer 




















Wood frog 


Rana sylvatica 


Thermal 


T 


Y(2,5) 






Y(S) 


Y(1,2) 




Freidenburg 






preference 
















and Skelly 






















(2004) 


Italian agile frog 


Rana latastei 


Larval mass 


T 


Y(2,5) 












Ficetola and 






Time to 




Y(2,5) 






Y(5) 


Y(1) 




Bernardi (2005) 






metamorphosis 






















Mass at 




N(2,5) 


















metamorphosis 


















Common 


Alytes 


Survival to 


P&T 




N(2,3) 










Richter-Boix 


midwife toad 


obstetricans 


metamorphosis 


(HY) 














et al. (2006) 






Time to 






Y(2,3) 




Y(5) 


Y(1) 










metamorphosis 






















Mass at 






Y(2,3) 




Y(5) 


Y(1) 










metamorphosis 


















Common parsley 


Pelodytes 


Survival to 






N(2,3) 












frog 


punctatus 


metamorphosis 






















Time to 






Y(2,3) 




Y(5) 


Y(1) 










metamorphosis 






















Mass at 






Y(2,3) 




Y(5) 


Y(1) 










metamorphosis 


















Common toad 


Bufo bufo 


Survival to 






N(2,3) 
















metamorphosis 






















Time to 






N(2,3) 
















metamorphosis 






















Mass at 






Y(2,3) 




Y(5) 


Yd) 










metamorphosis 


















Natterjack toad 


Bufo calamita 


Survival to 






N(2,3) 













metamorphosis 
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Common Name 



Species 



Trait Type 



Factor Genetic Plastic GxE 



Adapt 



Cause Time Reference 



Mediterranean 
tree frog 



Perez's frog 



Marsh frog 
Pool frog 
Edible frog 
Red back 
salamander 
Italian newt 
Italian 

crested newt 
Common frog 



Common toad 



Common toad 

Tokyo 

salamander 
Montane brown 

frog 

Forest green 
treefrog 
Palmate newt 



Hyla meridionalis 



Rana perezi 



Pacific treefrog Pseudacris regilla 



Common frog 

Common frog 

Hourglass 
treefrog 



Rana ridibunda 
Rana lessonae 
Rana esculenta 
Plethodon 
cinereus 
Triturus italicus 
Triturus carnifex 

Rana temporaria 



Bufo bufo 



Bufo bufo 

Hynobius 

tokyoensis 
Rana ornativentris 

Rhacophorus 

arboreus 
Lissotriton 

helveticus 



Rana temporaria 

Rana temporaria 

Dendropsophus 
ebraccatus 



Time to 
metamorphosis 
Mass at 

metamorphosis 
Survival to P&T 

metamorphosis (HY) 
Time to 

metamorphosis 
Mass at 

metamorphosis 
Survival to 

metamorphosis 
Time to 

metamorphosis 
Mass at 

metamorphosis 
Survival T 



Body length P&T 

Body length 

Body length 

Color phenotype T 

Time to hatching T 
Time to hatching 

Spawning date T 
Spawning date 
Spawning date 
Spawning date 
Migration to 
ponds 
Body condition T 
Survival 

Spawning date T 

Spawning date 

Spawning date 

Number of eggs T 
Oviposition 
period 

Hatching success 
Time to hatching 
Oviposition rate 
Hatching rate 
Spawning date T 



Melanism 
Melanism 
Stage at hatching 
Larval mortality 



Yt 



T(L) 

T 

P 

(HD) 



Y(5) 



N(2,5) 



N(2,3) 
N(2,3) 
N(2,3) 
Y(2,3) 
Y(2,3) 
Y(2,3) 
Y(2,3) 
N(2,3) 
Y(4,5) 



Y(4) 
Y(4) 
Y(4) 



Y(2) 
Y(2) 

Y(4) 
Y(4) 
Y(4) 
Y(4) 
Y(4) 

Y(4) 
Y(4) 
Y(4) 

Y(4) 

Y(4) 

Y(2,3) 
Y(2,3) 

Y(2,3) 
Y(2,3) 
N(2,3) 
N(2,3) 
Y 

Y(6) 

Y(2,3) 

Y(2,3) 

Y(2,3) 

Y(2,3) 



Y(5) 
Y(5) 
Y(5) 
Y(5) 

Y(5) 



Y(5) 

Y(5) 
Y(5) 

Y(5) 
Y(5) 
Y(5) 
Y(5) 
Y(5) 

N(B) 
N(5) 
Y(5) 

Y(5) 

Y(5) 



N(2,5) 



Yd) . 

Yd) 

Yd) 

Yd) . 
Y(1 ,2) FD 



Yd ,2) FD 

Y(1 ,2) FD 

Y(1 ,2) FD 

Y(1) FD 



Y(1) 
Yd) 

Yd ,2) 
Yd ,2) 
Yd ,2) 
Yd ,2) 
Yd ,2) 

Yd) 
Yd) 
Yd ,2) 



FD 
FD 
FD 
FD 
FD 

FD 
FD 
FD 



Yd ,2) FD 



Yd ,2) FD 



N(5) Y(1) 
Y(1) 

N(5) Y(1) 
Y(1) 



Y(4,5) Y(1 ,2) FD 



Y(5) Y(1 ,2) 

Y(1) 

N(5);Y(5) Y(1) 
N(5);Y(5) Y(1) 
Y(5) Y(1) 



Richter-Boix 
et al. (2006) 



Govindarajulu 

and Anholt 

(2006) 
Tryjanowski 

et al. (2006) 

Gibbs and 

Karraker (2006) 
D'Amen 

et al. (2007) 

Sparks 
et al. (2007) 



Reading (2007) 

Kusano and 
Inoue (2008) 



Galloy and 
Denoel (2010) 



Phillimore 

et al. (2010) 
Alho et al. 

(2010) 
Touchon and 

Warkentin 

(2010) 
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Common Name Species 



Trait Type 



Factor Genetic Plastic GxE Adapt Cause Time Reference 



Pool frog 



Rana lessonae 



Bullfrog 

Panamanian 

golden frog 
Hokkaido 

salamander 

Common frog 
Bullfrog 

Northern cricket 
frog 

Lowland leopard 
frog 

Dwarf 

salamander 
Marbled 

salamander 
Tiger salamander 

Ornate chorus 
frog 

Mole salamander 

Southern toad 

Eastern 
narrowmouth 
toad 

Spring peeper 

Southern leopard 
frog 

Eastern spadefoot 
toad 
Wood frog 
Spring peeper 

Northern leopard 
frog 



Rana catesbeiana 

Atelopus zeteki 

Hynobius 
retardatus 

Rana temporaria 

Rana catesbeiana 

Acris crepitans 

Rana yavapaiensis 



Eurycea 

quadridigitata 
Ambystoma 

opacum 
Ambystoma 

tigrinum 
Pseudacris ornata 

Ambystoma 
talpoideum 
Bufo terrestris 

Gastrophryne 
carolinensis 

Pseudacris 
crucifer 
Rana 

sphenocephala 
Scaphiopus 

holbrookii 
Rana sylvatica 
Pseudacris 

crucifer 
Rana pipiens 



Predation 

vulnerability 
Size at 

metamorphosis 
Survival 
Time to 

metamorphosis 
Mass at 

metamorphosis 
Size at 

metamorphosis 
Larval growth 

rate 
Body size 

Chytrid Bd 

infection 
Body size 
Time to 

metamorphosis 
Leg length 

Chytrid Bd 

infection 
Chytrid Bd 

infection 
Chytrid Bd 

infection 

Median arrival 
date 

Median arrival 
date 

Median arrival 
date 

Median arrival 
date 

Median arrival 
date 

Median arrival 
date 

Median arrival 
date 

Median arrival 
date 

Median arrival 
date 

Median arrival 
date 

Peak calling date 
Peak calling date 

Peak calling date 



Y(2,3) 



T N(2,5) N(2 

Y(2,5) Y(2 

N(2,5) Y(2 

N(2,5) Y(2 

Y(2,5) Y(2 

T(E) . Y(6i 

T . Y(6i 

T . Y(2 

Y(2 

T(L) Y(2,5) Y(2 

T Y(2 
Y(2 

T . Y(6] 



T&P 



Y(4 
Y(4 
Y(4 
Y(4' 
N(4 

N(4; 

N(4. 

N(4! 

N(4. 

N(4! 

Y(4 
Y(4 

Y(4) 



3) N(2,5) . 
3) Y(2,5) Y(5) 

3) N(2,5) . 

3) N(2,5) . 

3) Y(2,5) . 

Y(5) 

Y(5) 

Y(5) 

Y(5) 
Y(5) 
Y(5) 

Y(5) 
Y(5) 
Y(5) 
Y(5) 



Y(S) 
Y(5) 

Y(5) 



Yd) 

Y(1) 
Yd) 
Y(1) 
Yd) 

Y(1, 2) . 

Y(1, 2) FD 

Yd) 
Y(1) 

Yd, 2) . 

Yd, 2) . 

Y(1, 2) . 

Yd, 2) . 

Y(1, 2) FD 
Yd, 2) FD 
Y(1, 2) FD 
Yd, 2) FD 
FD 
FD 
FD 

FD 
FD 
FD 

Yd, 2) FD 
Yd, 2) FD 

Y(1, 2) FD 



Orizaola 
et al. (2010) 



Liu et al. 

(2010) 
Richards- 

Zawacki (2010) 
Michimae (2011) 



(Alho et al. 
2011) 

Chatfield and 
Richards- 
Zawacki (201 ' 

Forrest and 
Schlaepfer 
(2011) 

Todd et al. 
(2011) 



Walpole et al. 
(2012) 
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Common Name 


Species 


Trait Type 


Factor 


Genetic 


Plastic 


GxE 


Adapt 


Cause 


Time 


Reference 


Green frog 


Rana damitans 


Peak calling date 






N(4) 








FD 




Bullfrog 


Rana cateseiana 


Peak calling date 






N(4) 








FD 




American toad 


Bufo americanus 


Peak calling date 






N(4) 








FD 




Boreal chorus 


Pseudacris 


Mortality 


P&T 


N(2,5) 


N(2,3) 


N(2,5) 








Amburgey 


frog 


maculata 


Time to 

metamorphosis 
Size at 

metamorphosis 


(HY) 


Y(2,5) 
Y(2,5) 


N(2,3) 
N(2,3) 


N(2,5) 
N(2,5) 




Yd) 
Yd) 




et al. (2012) 


Natterjack toad 


Bufo calamita 


Breeding habitat 
selection 


T(L) 




Y(6) 




Y(5) 


Yd ,2) 




Rannap 
et al. (2012) 


Gunther's toadlet 


Pseudophryne 
guentheri 


Embryonic 
survival 


P 

(SM) 


Y(2) 


Y(2) 


Y(2,5) 


N(B) 


Yd) 




Eads 

et al. (2012) 






Time to hatching 




Y(2) 


Y(2) 


Y(2,5) 


N(5) 


Yd) 










Body size 




Y(2) 


Y(2) 


Y(2,5) 


N(5) 


Yd) 






Striped marsh 


Limnodynastes 


Age at hatching 


T 




Y(2,3) 




Y(5) 


Yd) 




Niehaus 


frog 


peronii 


Size at hatching 
Larval growth 
rate 
Larval 

developmental 
rate 

Age at hatching 
Size at 

metamorphosis 
Mass at 

metamorphosis 


TV 




Y(2,3) 
Y(2,3) 

Y(2,3) 

N(2,3) 
Y(2,3) 

Y(2,3) 




Y(5) 
Y(5) 


Yd) 
Yd) 

Yd) 

Yd) 
Yd) 




et al. (2012) 


Common frog 


Rana temporaria 


Activity level 


T(L) 


Y(2,5) 


Y(2) 


Y(2,5) 


Y(5) 


Yd ,2) 




Orizaola 






Body 




Y(2,5) 


Y(2) 


N(2,5) 




Yd ,2) 




et al. (2013) 






morphology 






















Larval period 




Y(2,5) 


Y(2) 


N(2,5) 


Y(5) 


Yd ,2) 










Mass at 




Y(2,5) 


Y(2) 


N(2,5) 




Yd ,2) 










metamorphosis 






















Growth rate 




Y(2,5) 


Y(2) 


N(2,5) 


Y(5) 


Yd ,2) 







A 'Y' indicates that evidence was found for genetic or plastic responses in traits or that adaptability or causality was investigated; 'N' indicates evi- 
dence was not found; '.' indicates that it was not investigated. Numbers next to a 'Y' or 'N' denote the method of investigation invoked. FACTOR 
(climate factor proxy): CD -climate data, E- elevation, GT- gestation temperature, HD - hydration, HY - hydroperiod, IT- incubation temperature, 
L - latitude, NH - nest habitat, NT - nest temperature, P - precipitation, SM - soil moisture, SST - sea surface temperature, T - temperature, TV - 
temperature variability; GENETIC categories: 1 - animal models, 2 - common garden studies, 3 - comparison to model predictions, 4 - experimental 
evolution, 5 - space for time substitution, 6 - molecular genetic approaches; PLASTIC categories: 1 - animal models, 2 - common garden studies, 3 
- experimental studies, 4 - fine-grained population responses, 5 - individual plasticity in nature, 6 - space for time substitution; GxE categories: 1 - 
animal models, 2 - common garden studies, 3 - comparison to model predictions, 4 - experimental evolution, 5 - space for time substitution; 
ADAPT categories: 1 - reciprocal transplants, 2 - phenotypic selection estimates, 3 - genotypic selection estimates, 4 - Q st f st comparison, 5 - com- 
mon sense; CAUSE categories: 1 - common sense, 2 - phenotype by environment interactions, 3 - experimental selection/evolution; for full descrip- 
tions of all categories see Merila and Hendry (this volume); TIME (time component included in data collection): RS - resurrection study, EX - field or 
greenhouse experiment through time, FD - field observations through time, MD - modeled through time. 
*Feature of temperature-dependent sex determination, 
f Other studies show color phenotype has a genetic basis. 



to warmer climates (Beebee 1995; Blaustein et al. 2001; 
Schwanz and Janzen 2008; Telemeco et al. 2009; Todd 
et al. 2011). In general, the phenological responses of 
amphibian breeding are among the largest observed - esti- 
mated to be two to four times stronger than those in other 
taxonomic groups (Parmesan 2007). Thus, the plastic 
response to climate change is well established and strong in 



amphibians. Generally, these rapid changes in phenology 
are assumed to occur too quickly to entail evolution. How- 
ever, the genetic component of this response has never been 
evaluated in amphibians or reptiles in common garden 
experiments. One study estimated no additive genetic vari- 
ance for breeding phenology in painted turtles, but moder- 
ate genetic variation emerged following a warm winter 
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Figure 1 Word clouds indicating most commonly analyzed climate- 
change-related traits for (A) amphibians and (B) reptiles. Traits are listed 
alphabetically from left to right with font size set to the number of 
times the trait was studied relative to others. 
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(V A = 0.13; evolvability CV A of 0.23; McGaugh et al. 
2010). Another study evaluated variation in temperature- 
dependent breeding dates along a latitudinal gradient and 
used observed differences in reaction norms to argue for 
the evolution of breeding phenology (Phillimore et al. 
2010). Although the differences in reaction norms are com- 
pelling, other uncontrolled environmental factors could 
explain these differences. Applying complicated statistics to 
observational data cannot replace the simple elegance of 
the common garden experiment. Some of these phenologi- 
cal changes over longer time periods might involve genetic 
changes, but we lack clear evidence to date. Introducing a 
population from another climate region into a common 
environment and evaluating differences in phenology 
between introduced and native populations under the same 
climate cues would indicate if these responses diverge 
genetically across landscapes. Determining if they evolve in 
response to climate change is more difficult, requiring 
researchers to identify candidate genes, parameterize 
genetic models, or experimentally estimate trait changes 
through time. 

Temperature-dependent sex determination presents 
another important case of plasticity with strong links to cli- 
mate (Janzen 1994; Mitchell and Janzen 2010). Warmer cli- 
mates generally bias sex ratios in reptiles toward females. 
Some populations could face demographic collapse if high 
temperatures preclude the development of one sex (Janzen 
1994). Yet, multiple studies suggest that reptiles can 
respond behaviorally through nesting site choice, depth, 
and breeding time (Doody et al. 2006; Doody 2009; Refsn- 
ider and Janzen 2012). These behaviors could mediate the 
effect of hotter temperatures on sex ratios. For instance, 
water dragons (Physignathus lesueurii) maintain similar 
incubation temperatures despite changing environmental 
temperatures by digging shallower nests and selecting more 
open nest sites in colder regions (Doody et al. 2006; Doody 



(B) 100 



c 

CD 
o 

i_ 
CD 




Adaptive plasticity 

Figure 2 (A) Percent of traits from weighted proportional meta-analy- 
sis (±95% confidence intervals) that show significant genetic, environ- 
ment, and genotype-by-environment contributions to responses to 
climate change and (B) percent of studies recording adaptive plasticity 
for amphibians (gray) and reptiles (white). Individual species and popu- 
lations within a given study are treated as separate estimates. Numbers 
within bars denote the total number of examples that analyzed a given 
phenotypic mechanism. 



2009). In contrast, changes in breeding time and nest depth 
did not buffer an observed 1.6°C warming period for the 
lizard Bassiana duperreyi. These lizards nested 4 weeks ear- 
lier and buried eggs 15 mm deeper, but eggs still warmed 
1.5°C, which threatened their normal sex ratios (Telemeco 
et al. 2009). Behaviors might moderate changes in sex 
ratios, but not always to the degree necessary to mitigate 
climate change. 

Different thermal regimes were associated with plastic 
changes in life-history traits such as clutch size, develop- 
ment rate, growth rate, and survival in amphibians (Voss 
1993; Galloy and Denoel 2010; Orizaola et al. 2010; Nie- 
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haus et al. 2011). Some of these changes were dramatic. 
For instance, half as many eggs were laid by palmate newts 
(Lissotriton helveticus) when raised at 22 vs 18°C (Galloy 
and Denoel 2010). Experimentally induced shorter hyd- 
roperiods altered survival, development, and growth in 
multiple European frogs and toads (Richter-Boix et al. 
2006). In contrast, simulated drying did not induce plastic 
changes in boreal chorus frogs (Pseudacris maculata; Am- 
burgey et al. 2012). 

Results suggest that plasticity will generally buffer some 
of the effects of climate change. However, not all traits are 
plastic, and sometimes plasticity might prove insufficient 
by itself to buffer climate change. 

Evidence for genetic responses to climate variation 

Our sample size was small and only one study explicitly 
evaluated genetic changes over time in response to climate 
change. However, all 17 studies that evaluated genetic dif- 
ferences found support for adaptive genetic differences in 
at least one trait. The estimated weighted proportion of 
genetic responses to climate variation in amphibians was 
0.65 (95% confidence interval: 0.46, 0.81), and the propor- 
tion in reptiles was 0.66 (0.43, 0.83; Fig. 2A). The un- 
weighted and weighted estimates did not differ 
substantially from each other, and random effects were not 
significant for amphibians and reptiles. 

In one amphibian example of a genetic response to cli- 
mate variation, wood frogs (Rana sylvatica) evolved differ- 
ent critical thermal maxima and thermal preferences 
between warm and cold ponds that had diverged in tem- 
perature over 36 years (Skelly and Freidenburg 2000; Frei- 
denburg and Skelly 2004). In another case, warmer climates 
led to increasing proportions of dark morphed red-backed 
salamanders (Plethodon cinereus; Gibbs and Karraker 
2006). Pool frog tadpoles (Rana lessonae) diverged in tem- 
perature-dependent developmental and growth rates across 
habitats that differed in temperature (Orizaola et al. 2010). 
For reptiles, snapping turtles (Chelydra serpentina) from 
different latitudes differed in their temperature-dependent 
sex determination patterns when raised in the laboratory 
(Ewert et al. 2005), and sex ratios differed among families 
for alligators and two turtle species (Rhen and Lang 1998). 
The overall magnitude of these changes was often not very 
large. For instance, the study on the evolution of critical 
thermal maxima in wood frog tadpoles found a difference 
of 0.4°C between populations inhabiting warm and those 
inhabiting cold ponds (Skelly and Freidenburg 2000). Simi- 
larly, the occurrence of striped red-backed salamander 
morphotypes decreased by only 6% over the last century 
(Gibbs and Karraker 2006). 

We only found six amphibian and three reptile estimates 
of genotype-by-environment interactions in response to 



climate variation. The reptile study (Rhen and Lang 1998) 
found evidence for genotype-by-environment interactions 
in two of the three species they examined. For amphibians, 
we had enough studies to apply a weighted meta-analysis. 
The estimated proportion of genotype-by-environment 
responses to climate variation was 0.48 (95% confidence 
interval: 0.13, 0.85; Fig. 2A). The unweighted estimate was 
similar (= 0.47), and random effects were not significant. 
One genotype-by-climate interaction occurred when 
R. sylvatica tadpoles varied in hatching depending both on 
if they came from warm or cold ponds and their exposure 
to cool or warm environmental conditions (Skelly and 
Freidenburg 2000). In R. lessonae tadpoles, larval growth 
rate and time to metamorphosis were accelerated the most 
in warm conditions for the population originating from 
the coldest environment, indicating countergradient varia- 
tion (Conover and Schultz 1995; Orizaola et al. 2010). 

Other research suggests the potential for genetic varia- 
tion to fuel the evolution of sex determination thresholds 
in reptiles in response to warmer climates. A review of rep- 
tiles with temperature-dependent sex determination 
reports heritabilities ranging from 0.26 to 0.82 (McGaugh 
and Janzen 2011). Morjan (2003) parameterized a model to 
predict that threshold temperatures might evolve more 
quickly than nest site choice, but not rapidly enough to 
counter predicted climate change. Morjan (2003) finds 
even slower rates of evolution of nest site choice owing to 
low heritability and high maternal contributions. Heritabil- 
ity of nest site choice in painted turtles (Chrysemys picta) 
ranged from 0.06 to 0.7 (McGaugh et al. 2010). Additive 
genetic variance increased in this trait following warm win- 
ters (McGaugh et al. 2010), which demonstrates that 
extreme future climates might reveal additive genetic varia- 
tion to fuel adaptation exactly when it is needed most. 
However, in other cases, limits to additive genetic variation 
might exist that would prevent evolution in response to 
more extreme climate variation. More generally, we high- 
light the many opportunities for researchers to perform 
novel experiments that test for the genetic basis of climate- 
related traits in amphibians and reptiles. 

On a cautionary note, few studies controlled for trans- 
generational effects, which could affect responses to climate 
change. Transgenerational effects occur when nongenetic 
contributions originate from previous generations (e.g., 
maternal effects). In addition, all our evidence for evolu- 
tion arose from genetic differences measured across popu- 
lations rather than through time. This space-for-time 
substitution is criticized as much as it is applied, but per- 
sists because of limited data. The problem is that past adap- 
tive divergence among populations need not indicate the 
potential for future adaptation, especially at the rates 
required by climate change. Spatial differences might have 
evolved over centuries, whereas climate change might 
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require evolution over decades. Also, strong natural selec- 
tion could produce demographic collapse before popula- 
tion-level fitness can recover (Gomulkiewicz and Holt 
1995; Sinervo et al. 2010). 

On a more positive note, the studies on R. sylvatica tad- 
poles revealed genetic differences between ponds that had 
undergone known temporal changes in temperature over 
just a few decades (Skelly and Freidenburg 2000; Freiden- 
burg and Skelly 2004), suggesting that space can substitute 
for time when we have information on recent environmen- 
tal changes across the landscape. At least some evolutionary 
rates might be sufficient to meet the demands of the chang- 
ing climate. 

Are traits adaptive? 

We next assessed if trait differentiation was adaptive or 
maladaptive based on authors' conclusions and logic. In 
the latter case, we assumed a trait was maladaptive if it 
reduced fecundity or survival. Only a few studies examin- 
ing genetic variation also provided information that could 
be used to assess if changes were adaptive. Of those studies, 
all six amphibian studies and the one reptile study indi- 
cated an adaptive origin for genetic changes. 

Climate change factors generally induced adaptive plas- 
ticity (Fig. 2B). The weighted proportion of plastic 
responses to climate variation in amphibians that was con- 
sidered adaptive was 0.67 (95% confidence interval: 0.51, 
0.79). The random effect of species was not significant 
(P = 0.813), but study and population significantly 
explained additional variation (y\ = 26.1, 8.2, respectively; 
P < 0.005). The unweighted model indicated the same pro- 
portion. The proportion of plastic responses to climate var- 
iation in reptiles that was considered adaptive was 0.62 
(95% confidence intervals: 0.36, 0.82). The random effects 
of species, study, and population did not significantly 
explain additional variation. The unweighted model indi- 
cated the same proportion. For both amphibians and rep- 
tiles, allochronic studies did not show significantly 
different proportions of adaptive variation when compared 
against synchronic studies (P > 0.6). 

Although most observed plasticity was assumed to be 
adaptive, we also detected substantial maladaptive plas- 
ticity (Fig. 2B). Hotter temperatures enhanced learning 
performance in the lizard B. duperreyi (Amiel and Shine 
2012) and accelerated development and increased 
juvenile body condition in lacertid lizards (Monasterio 
et al. 2011). However, higher incubation temperatures 
produced lower hatching success, slower righting 
responses, poorer swim performance, and skewed sex 
ratios in other reptiles (Janzen 1994; Wapstra et al. 
2009; Mickelson and Downie 2010; Micheli-Campbell 
et al. 2011; Bell et al. 2013). Despite the near- 



universality of plastic responses to climate variation, 
some plasticity will ultimately prove maladaptive under 
altered climates. Understanding the fitness consequences 
of climate-induced plasticity requires researchers to 
delineate plastic responses in experiments, create artificial 
phenotypes independent of genotype (Waddington 
1953), and then track their long-term fitness in experi- 
mental or natural situations. 

Which selective force? 

We also evaluated the degree to which studies evaluated the 
causal linkage between climate change and each phenotypic 
response. Although multiple methods exist (Merila and 
Hendry 2014), our studies assessed causality using com- 
mon sense, phenotype-environment correlations, or exper- 
imental manipulations. Climate was suggested to be 
involved in phenotypic changes in all the traits analyzed 
(Table 1). However, these conclusions were often based on 
common sense rather than rigorous tests. 

Future directions 

Given limited information on adaptive responses to climate 
change, we argue that we need to intensify our efforts to 
explore how plasticity and genetic adaptation will attenuate 
the future impacts of climate change. These mitigating fac- 
tors might not provide a panacea, but we should estimate 
their contributions and include them in any rigorous 
assessment of climate change effects on species. Most of the 
examples we found required substituting space for time 
and thus constituted weaker inferences. Quantitative 
genetic models might be used to predict future phenotypes, 
but they assume that we can accurately estimate the herit- 
abilities of traits and strength of selection in the wild (En- 
dler 1986), which is often not the case (Hoffmann and 
Merila 1999; Charmantier and Garant 2005). Genetic cor- 
relations can constrain evolutionary responses to climate 
change (Etterson and Shaw 2001). Moreover, attempts to 
predict evolution in the wild have frequently failed. These 
shortcomings are due to (i) poor estimates of heritability 
or selection in the wild, (ii) the ability of selection to act on 
nonheritable phenotypic variation, (iii) the fact that traits 
are genetically correlated, (iv) rapidly shifting environ- 
ments, or (v) because of a simple lack of statistical power 
(Merila et al. 2001). All of these reasons cast doubt on our 
ability to extrapolate current evolutionary responses far 
into the future. 

Other options include using analogous changes in land- 
scapes as natural experiments (Skelly and Freidenburg 
2000) or performing common garden experiments 
through time, where relevant phenotypic responses are 
evaluated in the same populations and using the same 
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methods across long enough periods to capture climate 
change. These rigorous experiments would reveal evolu- 
tion in real time, highlight if evolution is keeping up with 
changing climate, and provide insights into the limits of 
additive genetic variation to fuel evolutionary responses in 
the near future. Experimental evolution in the wild pro- 
vides another compelling option. Researchers could design 
experiments that manipulate climate factors directly and 
evaluate evolutionary changes in focal populations or 
transplant populations into selection regimes that mimic 
future climates and then track their evolution. However, 
the lack of allochronic studies in amphibians and reptiles 
is clearly related to the fact that these species often have 
long generation times that prevent such comparisons 
being made over practical time spans. Lastly, a candidate 
genes approach offers many advantages if sufficient 
genomic and trait-gene mapping is available or can be 
accomplished. 

Conclusions 

Many amphibians and reptiles have already undergone 
climate-related extinctions or declines (Gibbons et al. 2000; 
Parmesan and Yohe 2003; Pounds et al. 2006; Sinervo et al. 
2010), and others are approaching the edge of their thermal 
tolerance (Huey et al. 2009). Models predict future declines 
in amphibians and reptiles owing to climate change (Tho- 
mas et al. 2004; Malcolm et al. 2006). However, these correl- 
ative models rarely account for evolutionary change. Strong 
plasticity or high genetic variance might prove sufficient to 
generate phenotypic changes that match the demands of a 
new climate for some species. Alternatively, relatively long 
generations and limited genetic variation could reduce their 
capacity for adaptive change after plastic reserves have been 
exhausted (Morjan 2003; Sinervo et al. 2010). 

Here, we show that both adaptation and plasticity could 
play important roles in mitigating climate change for 
amphibians and reptiles, but data is still limited, especially 
for genetic responses. Our literature review suggests that 
plasticity might often mediate changes to climate in both 
taxonomic groups. A long history of research has focused 
on plasticity in these species, and we show here that chang- 
ing climates might often induce adaptive plasticity. 
Although evidence for genetic adaptations to climate was 
limited to a handful of mostly synchronic studies, all of 
these studies found evidence for climate adaptation in at 
least one trait. Similar assessments of other taxonomic 
groups in this issue provide evidence for widespread plastic 
and adaptive responses to climate variation and concur- 
rently highlight the same lack of data for evolutionary 
responses (Franks et al. 2014; Reusch 2014; Schilthuizen 
and Kellermann 2014). If future studies confirm these 
results, then plasticity and adaptive capacity could dampen 
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some of the doom and gloom associated with research on 
biotic responses to climate change. 
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